Interplay between disorder and intersubband collective excitations in the 

two-dimensional electron gas 
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Intersubband absorption in modulation-doped quantum wells is usually appropriately described 
as a collective excitation of the confined two-dimensional electron gas. At sufliciently low electron 
density and low temperatures, however, the in-plane disorder potential is able to damp the collective 
modes by mixing the intersubband charge-density excitation with single-particle localized modes. 
Here we show experimental evidence of this transition. The results are analyzed within the frame- 
work of the density functional theory and highlight the impact of the interplay between disorder 
and the collective response of the two-dimensional electron gas in semiconductor heterostructures. 
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Intersubband (IS) excitations of two-dimensional elec- 
tron gases (2DEGs) confined in semiconductor het- 
erostructures are among the -most important probes of 
electron-electro]ji|rinteractionaJI and represent an active 
area of researchBu Further motivation for the experimen- 
tal and theoretical study of IS transitions stems from the 
growing number of opto-electronic-devices whose opera- 
tion is based on these transitions.alj In the last decade 
several authors have demonstrated that -particularly in 
the long-wavelength limit- IS charge-density modes (the 
ones probed in optical absorption experiments) are sig- 
nificantly shifted from the single-particle transition ener- 
gies by both direct {depolarization shift) and exchange- 
correlation J eEpifonic coupling) terms of the Coulomb 
interaction.cl Lj The influence of these dynamic contri- 
butions to the elementary excitation spectrum of 2DEGs 
has been studied experimentally and theoretically mod- 
eled. Only in recent times, however, attention has been 
given to the ipterBlay between disorder and many-body 
contributions.li2rii3 This interplay raises several funda- 
mental questions on the role of dephasing and scattering 
in determining the collective response of electrons and 
has direct implications on the optimization of IS-based 
devices. 

Metzner et a/.0 and Ullrich et al^ addressed the influ- 
ence on the IS linewidth of the various scattering mech- 
anisms and developed microscopic theories able to treat 
many-body effects and disorder on equal footing. One 
of the main results of these models is the precise de- 
scription of how the in-plane disorder destroys the co- 
herence of the collective excitations and modifies the IS 
linewidth and peak energy. In particular, Metzner et al. 
emphasized that a transition from a regime where the 
intersubband absorption is rather broad and originates 
from randomly distributed localized single-particle exci- 
tations to a regime dominated by a sharp and blue-shifted 
collective mode of the 2DEG can occur by varying the 
electron density. This effect was explained as a mutual 
phase adaptation of the localized oscillators interacting 



via long-range Coulomb forces. Recently Yakimov et al. 
reported experimental evidence of such a transition in 
the case_pf holes in a self- aggregated Si-Ge quantum dot 
matrix.Ej 

In this article we report a study of mid-infrared IS 
excitations in 2DEGs confined in modulation-doped nar- 
row GaAs/Alo.sGao.yAs quantum wells (QWs) under the 
application of a perpendicular electric field. The latter 
allows us to vary the carrier sheet density in the range 
10^-10^^ cm~^ and explore regions where the relative 
role of disorder and many-body effects are expected to 
be markedly different without reaching the extreme lo- 
calization limit studied by Yakimov et al. in Ref. We 
present quantitative analysis of IS transitions based on a 
model that includes both static and dynamic many-body 
effects. Our measurements reveal a substantial softening 
of the IS excitation energy accompanied by an abrupt 
increase of the intersubband absorption linewidth. This 
behavior is observed at low temperatures and at elec- 
tron densities ~10^^ cm~^. We shall argue that at these 
densities a transition occurs from a regime dominated by 
the collective response of free electrons to another driven 
by the influence of localized single-particle intersubband 
modes on this response. This is analogous to the Landau 
damping mechanism occurrins-at finite in-plane wavevec- 
tor in a disorder-free 2DEGJa't3 

The two samples used for this study are modulation- 
doped single GaAs/Alo.sGao.yAs QWs (8.7-nm- and 7.5- 
nm-thick, respectively). Doping in the Alo.3Gao.7As 
layers was offset by about 25 nm from the wells, and 
the growth was performed by molecular beam epitaxy. 
At equilibrium samples have free-electron concentrations 
close to 10^^ cni~^ and low-temperature mobility above 
10^ cm^/Vs. Well thickness was chosen in order to avoid 
population of the second subband and to maximize the 
impact of fluctuations due to interface roughness, while 
preserving high values of electron mobility. Measured 
low-temperature IS transitions peak at 112.5 meV (8.7- 
nm-thick QW) and 130.5 meV (7.5-nm-thick QW). In or- 
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der to vary continuously the electron density in the sam- 
ples, a metallic gate was evaporated on the surface. In 
this way intersubband absorption spectra were measured 
in the carrier-concentration range 2x 10^-10^^ cm~^ 
and in the temperature range 5-300 K. The 2DEG den- 
sity was carefully determined by Shubnikov-de Haas and 
Hall measurements as a function of gate voltage Vg. For 
optical measurements, we fabricated 45°-edge multipass 
waveguides. Absorption spectra were acquired with an 
infrared Fourier transform interferometer in normal or 
step-scan modes and using infrared light polarized per- 
pendicularly to the QW layers and propagating under the 
gate. In the step-scan case, the bias voltage was square- 
wave modulated between -2V and the desired voltage, at 
a frequency between 1 and 20 kHz. In both cases, back- 
ground spectra were collected at Vg = — 2 V, at this bias 
the QWs were fully depleted. 

Figure |l] shows representative absorption spectra for 
the 8.7-nm- (left panel) and 7.5-nm-QW structures (right 
panel) out of a large set taken at many gate-bias val- 
ues and different temperatures. Dashed lines refer to 
electron density n « 8.4 x 10^^ cm~^ {Vg = 0.1 V and 
0.4 V in the 8.7 and 7.5 nm QWs, respectively), sohd 
lines to n « 4x10^° cm"^ [Vg = -1.2 V, left panel) and 
n « 1x10^° cm"2 (y^ = i right panel). Data 
correspond to transitions between the first two subbands 
in the well, as depicted in the inset of Fig. |l|. At the 
higher electron densities, spectra have a lorentzian line- 
shape with a full width at half maximum (FWHM) of 
3.5 and 4.5 meV. Surprisingly, the measured peak-energy 
shift at the different biases in Fig. |l|is very small, despite 
the very different electric fields (up to 8x10''^ V/cm at the 
lower carrier densities) and carrier concentrations. In- 
deed, a single-particle description for the IS energy that 
includes only the static Hartree term (and therefore takes 
into account exclusively the appropriate band-bending of 
the biased heterostructures) does lead to the prediction 
of a significant energy shift (see Fig. ^, dashed line for the 
case of the 8.7 nm QW). The observed weak dependence 
of IS transition energy on electron density is the result 
of a cancellation effect between Stark shift and dynamic 
many-body contributions. El This is clearly shown in Fig. ^ 
where measured IS peak energies for the 8.7 nm structure 
(dots) are plotted as a function of carrier sheet density to- 
gether with results of a theoretical model (solid line) that 
includes many-body contributions within the adiabatic 
local density approximation (LDA) of the density func- 
tional theory, and therefore explicitly takes into account 
the effects of depolarization and excitonic shifts. The IS 
peak energy was calculated by evaluating the maximum 
of: 
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where dzzi^) is the zz component of the frequency de- 
pendent conductivity tensor, and G12 is proportional to 
the response function of a non-interacting 2DEG includ- 
ing the specific energy dispersions of the two individual 



subbands considered. These are evaluated extending the 
calculations reported in Ref. |l^ and considering the band 
bending in the entire heterostructure (see also Ref. |^) ; in 
all the theory differentojon-parabolic dispersions in every 
layer were considered. Il3 In (|l]) the depolarization contri- 
bution is linked to 



e{z',uj) 



£12 ' 



(2) 



where e{z,Lu) is the layer-dependent dielectric function, 
whose frequency dependence is estimated like in Ref. |^, 
^1^2 are the wave-functions along the growth direction 
and £12 is the difference between the minima of the two 
subbands, as calculated from a Schrodinger type equa- 
tion derived in the LDA formulation of density functional 
theory. The excitonic contribution is linked to 
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where v^cdN) is the exchange-correlation potentiallij as 
a function of the position-dependent three-dimensional 
electron density N determined self-consistently. 

The inclusion of all these effects was necessary to ob- 
tain the successful description of the experimental re- 
sults in the broad range of carrier density that is shown 
in Fig. 1^. In particular, the accuracy of our model in 
determining the 2DEG density is shown in the inset 
of Fig. |2[ where the carrier density estimated by Hall 
measurements (short-dash line), Shubnikov-de Haas os- 
cillations (empty circles) , and area under the absorbance 
peak (with oscillator strength calculated as in Ref. 0) IS 
reported as a function of Vg together with the theoretical 
dependence as derived within LDA (solid line) . Figure ^ 
shows also an excellent agreement between experimental 
and theoretical results for IS transition energies in a wide 
range of electron density (similar results were obtained 
for the second sample). This agreement however is not 
present at low densities where the theory is not able to 
reproduce the softening of the IS transition energy. This 
behavior signals a change in the nature of the intersub- 
band transition and originates from an interplay between 
many-body effects and disorder. This is the focus of this 
paper and in what follows we shall discuss the manifes- 
tations of this transition and its physical meaning. 

In order to analyze the low-electron density behavior 
of IS transitions, it is useful to refer to Fig. || where 
the deviation of the experimental data from the calcu- 
lated values is reported for both samples (open squares, 
right scale) at T = 5 K. At high density this deviation 
is «0 meV for both samples, but for n < 3x 10^^ cm^^ 
7 nm sample, upper panel) and n < 4 x 10^^ cm"'^ 
(7.5 nm sample, lower panel) we observe a notable drop 
in energy of about 2 meV and 3.5 meV, respectively. 
This evolution was monitored down to very low densities 
(1.4x10^" cm~^ and 2x10^ cm~^, respectively). The pro- 
nounced decrease of the intersubband absorption spectral 
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weight hindered measurements at even lower densities. 
Remarkably, the softening of the IS transition energies 
occurs together with a rather sudden increase of the cor- 
responding FWHM values. This is also shown in Fig. ^ 
(solid circles, left scale) where experimental FWHMs are 
plotted as a function of electron density. In the high den- 
sity regime the FWHM stays close to a quite small value 
(around 4 meV) . This type of behavior is consistent with 
calculations that include many-body contributions. In 
fact it has been shown, and confirmed by our calcula- 
tion, rtiiat these terms bring to its intrinsic homogeneous 
valuella the absorption-peak linewidth drastically reduc- 
ing it from the values calculated taking into account the 
different and non-pacaholic dispersions of the subbands, 
that increase with nJlj^E^ Both the energy position and 
FWHM trend, therefore, indicate that in the high-density 
regions absorption originates from purely collective inter- 
subband excitation (charge density excitation - CDE). 

The observed behavior at the lower densities signals 
that a change takes place in the collective response of 
the 2DEG. This transition occurs when the Coulomb in- 
teraction between electrons is not strong enough com- 
pared to disorder and cannot concentrate the whole os- 
cillator strength in a narrow CDE modeU In this regime 
the IS transition incorporates contributions of weakly- 
interacting randomly localized excitations, and this leads 
to the observed peak broadening and its shift towards 
lower energies. This evolution is schematically illustrated 
in Fig. ^, where the trends in the CDE energy (black line) 
and in the single particle excitation (SPE) continuum 
(darker gray areas) are shown as a function of density. 
The light gray area represents the intrinsic CDE broaden- 
ing. The darkest gray area includes the non-parabolicity 
contributions. The dominant source of broadening for 
single-particle excitations (SPEs), however, is inhomoge- 
neous in nature and gives rise to the overall SPE broaden- 
ing shown in Fig. |4[ The observed transition takes place 
when the broadened CDE and SPE excitations merge. 
At this density new relaxation channels of the intersub- 
band CDE into localized single-particle modes open, a 
phenomenon similar to corvYpjtional Landau damping at 
finite in-plane wavevectors.&EJ At even lower densities in- 
tersubband spectra are characterized by localized SPEs, 
This second regime was explored by Yakimov et al.tJ 
thanks to the large fluctuations provided by the self- 
aggregated quantum dots in the samples they studied. 
This overall behavior was theoretically discussed in Ref. 
|lT| . It can be noted that the transition region here dis- 
cussed occurs at higher density for the narrower well, 
where indeed the electron states are more sensitive to 
well-width fluctuations. Furthermore, we observed that 
these changes in energy and linewidth at low densities 
display a marked temperature dependence and were not 
detected at temperatures larger than 60 K (data not 
shown). This suggests an intriguing thermally-induced 
modification of the collective response of the system. It 
must be noted in fact that, contrary to this behavior, the 
many-body contributions to the peak position persisted 



at all temperatures studied. 

However, in order to further corroborate our inter- 
pretation, we developed a simple model based on the 
Drude approach for the response function of a QW char- 
acterized by a gaussian-broadened single-particle inter- 
subband transition with a phenomenological inhomoge- 
neous broadeniug, superimposed to a costant homoge- 
neous linewidth.tj Results for the deviation of the peak 
position from the depolarization-shifted energy and for 
the FWHM, as obtained by a lorentzian fit of the calcu- 
lated response function, arc in agreement with the data 
reported in Fig. ^. Within this simplified model, the 
transition between the purely collective excitation and 
the intermediate regime discussed above occurs when the 
depolarization shift is comparable to the inhomogeneous 
broadening consistently with our data. Also, the soft- 
ening of the transition energy occurs because the CDE 
peak, located in the high-energy end of the broadened 
single-particle excitation spectrum, possesses all the os- 
cillator strength at high density. 

In conclusion, we have experimentally and theoret- 
ically studied intersubband transitions in modulation- 
doped narrow GaAs quantum wells in a wide range of 
electron densities. We have reported evidence of the 
damping of the CDE due to coupling with single-particle 
localized modes as the density is lowered. This transition 
is a direct manifestation of the influence of the interplay 
between disorder and many-body effects in determining 
the nature of IS excitations of the 2DEG in semiconduc- 
tor heterostructures. 

We thank A. Tredicucci, G. Vignale and B. Vinter for 
useful discussion. The work at Scuola Normale was par- 
tially supported by MURST. 
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FIG. 2. Low-temperature intersubband absorption peak 
energies as a function of carrier sheet density for the sam- 
ple with 8.7 nm well width; experimental data are shown as 
dots with error bars. Solid line: calculated peak energy in- 
cluding many-body effects within the adiabatic local density 
approximation of the density functional theory. Dashed line 
is the Hartree approximation. Inset: electron sheet density as 
a function of applied bias as deduced from Hall (dotted line), 
Shubnikov-de Haas (open circles), and optical measurements 
(squares). Solid line is the result of the theoretical model. 
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FIG. 1. Low-temperature quantum well (QW) intersub- 
band absorption spectra at high (dashed line) and low (solid 
line) electron density for the two structures studied. Mea- 
surements were performed in a 2.6-mm-long and 370-/im-high 
optical waveguide. Dashed lines correspond to electron 
density n ~ 8.4 x 10^^ cm~^. Solid lines correspond to 
n « 4x10^° cm"^ (left panel), and to n ^ lxW^° cm"^ 
(right panel). Inset: Calculated band diagram for the 8.7 nm 
QW with no applied bias. Horizontal solid lines correspond to 
the energy of the bottom of the first two subbands; horizontal 
dashed line is the chemical potential. 
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FIG. 3. Full width at half maximum (FWHM, sohd cir- 
cles) of the intersubband absorption peak and deviation of 
measured peak position from theoretical values {AE, open 
squares) as a function of electron density. 
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FIG. 4. Schematic representation of charge-density excita- 
tion (CDE) energy (solid line) and intersubband single par- 
ticle excitation (SPE) continuum (gray and dark gray areas) 
versus electron sheet density. The light gray area represents 
the intrinsic CDE broadening. The darker gray area shows the 
broadening originating from subband dispersions, the gray re- 
gion the inhomogeneous broadening. 
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